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Abstract. We report on the fabrication of hybrid molecular devices employing multi¬ 
layer graphene (MLG) flakes which are patterned with a constriction using a helium 
ion microscope (HIM) or an oxygen plasma etch. The patterning step allows for the 
localization of a few-nanometer gap, created by electroburning, that can host single 
molecules or molecular ensembles. By controlling the width of the sculpted constriction, 
we regulate the critical power at which the electroburning process begins. We estimate 
the flake temperature given the critical power and find that at low powers it is 
possible to electroburn MLG with superconducting contacts in close proximity. Finally, 
we demonstrate the fabrication of hybrid devices with superconducting contacts and 
anthracene-functionalized copper curcuminoid molecules. This method is extendable 
to spintronic devices with ferromagnetic contacts and a first step towards molecular 
integrated circuits. 
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In recent years, molecular electronics has greatly benefited from the fabrication of 
nanometer spaced gold electrodes in the form of electromigrated and mechanical break 
junctions HD- While these advances have led to a wealth of new physics in molecular 
devices, there are improvements to be made and a continuing focus in the field is toward 
making robust and reliable contact to single molecules |8|. One recent direction is the 
use of graphene or multi-layer graphene (MLG) as an electrode and contacting molecules 
by covalent bonding or n — i r stacking of aromatic rings [9-12]. This method resulted in 
observation of the gating of a single molecule device at room temperature [13]. While gold 


atoms are mobile at room temperature [14,15 , MLG can be used as a stable intermediary 
material to contact molecules. In addition, two-dimensional MLG bridges better the size 
discrepancy between single molecules and bulky 3D electrodes where screening of the gate 
can reduce coupling to the molecule. Using MLG as an intermediary material also allows 
the option of contacting superconducting or ferromagnetic metals without changing the 
anchoring chemistry of the molecules. This opens the door to the fabrication of hybrid 
devices to study superconductivity and spintronics in single molecules. Up to now, the 
added complexity in fabricating hybrid superconducting-molecular devices has kept their 
exploration to a minimum [l6|. 

In this paper, we present further steps in the development of MLG electrodes with 
gold contacts and improve yield and stability through the localization of a constriction 
in an exfoliated flake. A predefined constriction using a helium ion microscope (HIM) 
to mill the flake or an oxygen plasma to etch the flake allows the flexibility to choose 
the location of a few-nm gap created by a process of electroburning. In comparison with 
as-exfoliated devices (without a constriction), MLG break junctions (with a constriction) 
show a 15% increase in yield due to a reduction in abrupt breaks and reconnected junctions 
after electroburning [l7|. The key to using MLG as opposed to graphene or even few- 
layer graphene is two-fold: the conductance of MLG is independent of back-gate voltage 
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making electron-transport measurements solely a study of the molecular characteristics of 


the device and not of the MLG electrodes, and we have found, as detailed in Ref. 17 


that the highest yield for nanogap formation is in devices having a resistance of 500- 
lkfl corresponding to devices with a MLG flake of 10 nm thickness with a variation 
of approximately 5 nm. Furthermore, we explore the possibility of fabricating hybrid 
superconducting-molecular devices. The constriction results in a significant reduction in 
the power required to start the burning process making it possible to fabricate MLG break 
junctions with superconducting contacts. We use the superconducting alloy, molybdenum- 
rhenium (MoRe 60/40) which has a higher melting temperature (Tm = 1150 °C) as 
compared with commonly used superconductors such as aluminum ( Tm = 660 °C) [l8j. 
Out of 26 fabricated MLG break junctions with MoRe contacts, 16 junctions resulted in 
tunnel gaps with an average low bias (Vj, = 10 mV) resistance of 17 ± 11 GO. 

A brief summary of the fabrication and electroburning process for MLG electrodes is 
given here which is detailed in Refs. 13,171. MLG devices are fabricated starting with 
heavily-doped Si wafers with a 285 nm-thick Si02 him. The Si substrate is used as a gate 
electrode for three-terminal devices. MLG hakes are deposited by mechanical exfoliation 
onto the Si 02 and hakes with a thickness of approximately 10 nm are chosen by their color 
contrast under an optical microscope. Typical hakes have a deep purple color under white 
light illumination and a red channel contrast of about 16% (see supplement information). 
Gold contacts, or MoRe contacts in the case of hybrid devices, are created by patterning 
with e-beam lithography and thin-hlm metal deposition. 

Electroburning is performed at room temperature and in air. A bias voltage up to 10 
V, corresponding to an electric held of more than 10 6 V/m for typical devices, is ramped 
up across the hake. At large current densities, removal of carbon atoms in the MLG hake 
occurs due to the interaction with oxygen triggered by the high temperatures of Joule self 
heating I3jl7|. Feedback controlled software is used to control the voltage across the hake. 
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Figure 1. (a) AFM image of a MLG device before electroburning. The saturated white 
regions at the top and bottom are gold electrodes, (b) AFM image of the same MLG 
flake after electroburning. The gap is seen near the top electrode indicated by a white 
arrow, (c) Critical power required to start the electroburning vs. the width of the MLG 
flake. Red circles show the critical powers for as-exfoliated devices. Blue circles show 
the critical power for pre-patterned break junction devices using the HIM. The dashed 
line shows the calculated critical power as a function of width (see text). Inset shows the 
calculated temperature profile for the widest flake of 4.5 p,m. 



When the conductance drops by 10% the voltage is swept down at a rate of 0.1 V//is to 
arrest the burning process and prevent the formation of a large gap. This cycle continues 
until the low bias (14 = 10 mV) resistance is greater than 10 Gf2, independent of sample 
geometry. This resistance ensures the electroburning of any last remaining nanometer¬ 
sized graphene islands which have been shown to act as large addition energy quantum 


dots 19 . By controlling the voltage in this way, a gap is burned across the entire flake 


with a width of a few nanometers at its closest point 20 


In Fig. 1 we show the electroburning results of as-exfoliated MLG flakes. Fig. 
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1(a) shows an atomic force microscopy (AFM) image of a few-micron-wide MLG flake 
before electroburning. Without localization of a constriction, the electroburning process 
is believed to most frequently start along the edge of the flake at an unsaturated carbon 
atom site where a lower displacement energy (7 eV) is required to remove an atom from 
the lattice as compared to removing a saturated atom from the flake center (30 eV) |21,221. 
In Fig. 1(b) the resulting few-nm gap is shown near the top of the electrode (white arrow). 
Electroburning of as-exfoliated flakes often results in gaps with unpredictable locations. 
Here, the close proximity of the gap to the gold electrode could result in screening of an 
applied gate field. Ideal devices should have gaps located near the center of the flake where 
the gate coupling to the molecule is unaffected. 

In Fig. 1(c), the critical power required to start the electroburning process for the 
flake shown in Fig. 1(a) is plotted along with several other flakes (red circles). Critical 
powers reach hundreds of rnW for the widest flakes and decrease with the flake width. 
Assuming comparable flake width (W) and length (L), the critical power as a function 
of width can be estimated by P c = I 2 R = ( j c Wt) 2 R s L/W ~ ( j c t) 2 R s W 2 , where we take 
j c = 5.8 x 10' A/cm 2 as the average critical current density in these devices 13 and 


R s = 100 fl/sq as the sheet resistance for MLG 23 . Using these values and an average 
flake thickness of t = 10 nm, we plot this critical power relation as a dashed line in 
Fig. 1(c) where the scatter in the data is mostly due to differences in flake geometry and 
thickness. 

The electroburning is triggered by the high temperatures from Joule heating of the 
MLG flake. It is interesting to estimate the temperatures reached given the critical 
power of the electroburning process for as-exfoliated MLG devices. An estimate of the 
temperature distribution along the flake (in the direction of current) can be made using 
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the ID heat equation 24 


4 ,€ 


+ p' x ~ d{T — T 0 ) — 0, 


( 1 ) 


where A = Wt is the cross-sectional area of the flake, k is the thermal conductivity of 
the MLG flake, p' x is the Joule heating rate in watts per unit length, and g is the thermal 
conductance to the substrate per unit length and is calculated given the thermal resistance 
of the Si0 2 and Si substrate beneath the MLG flake (g ~ l/[L(R ox + Rsi)}, where L is 
the length of the flake (in the direction of current), R ox = t ox /{k ox WL ), where k ox = 1.4 

. Heat convection 


WK -1 m -1 and R Si = 1/{2k si (WL) l / 2 ) } where k Si = 50 WK’ 1 


m 


25 


to air is neglected as the heat convection coefficient for air (10 — 100 Wm _2 K _1 ) is four 
orders of magnitude smaller than the coefficient for heat conduction to the substrate 
(calculated below). 

An analytic solution to the ID heat equation is arrived at by assuming that k is 
constant and that there is a uniform heat generation along the flake given by p' x ~ I 2 R/L. 
The temperature in the flake at position x is then given by: 


rpf \ nr , Px i COsh{x L h ) 

T{x) = T 0 + — 1- ut /or ' 

g \ cosh(L 2L h ) 


( 2 ) 


Here, T 0 = 300 K, is the temperature of the electrodes and Lh = ( kA/g) 1 ^ 2 is the 
characteristic healing length along the flake and is a measure of lateral heat dissipation 
in the flake. Using Eqn. [2j we estimate the critical temperature distribution for the 
widest as-exfoliated flake in Fig. 1(c) (W = 4.5 //m, L — 2 /irn) with a critical power of 
nearly 190 rnW. We take the thermal conductivity of MLG (k = 1000 WK^ivT 1 ) which 


is found for flakes of 8 layers 27 . From the length and width of the flake, we calculate 
the thermal resistance of the Si substrate R$i ~ 3 x 10 3 K/W, the thermal resistance of 
the oxide layer R ox w 2x 10 4 K/W, the thermal conductance per unit length g & 19 
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WK 1 m x , and the characteristic healing length for this device, Lh = ( kA/g ) 1//2 = 1.5 
//ni. The healing length is larger than estimates for graphene (Lh ~ 0.2 /irn on Si02) 


due to the thicker cross sectional area 25 . The temperature distribution is plotted in the 
inset of Fig. 1(c). The center of the flake is expected to reach temperatures of ~ 1200 K. 
Typical device dimensions are larger but comparable to the healing length making the ID 


diffusion equation a sufficient estimate of the temperature distribution in the flake 25 


A more accurate calculation can be made using finite element methods with COMSOL 
MULTIPHYSICS v4.3b to account for the small lateral heat dissipation (see supplemental 
information). The high flake temperatures prove detrimental to the fabrication of MLG 
electrodes with MoRe contacts, demonstrated below. These temperatures can be reduced 
by pre-patterning the flake to make MLG break junctions with decreased widths. 

The electroburning of MLG break junctions is first demonstrated with gold contacts 
in Fig. 2. Patterning of the MLG devices before electroburning is achieved using a HIM 
or by oxygen plasma etch. For oxygen plasma etching, we use a pressure of 50 /iBar with 
a flow of 25 seem O 2 at 20 W power for 30 secs. The HIM is well suited for high resolution 
nanofabrication and unlike oxygen plasma etching, HIM patterning does not require an 


e-beam lithography step which makes the fabrication process cleaner and faster 28 . In 
comparison with gallium ion beam milling, HIM patterning allows higher resolution and 


results in less damage to the Si/Si02 substrate 29 . Patterning of graphene using a HIM 


has been achieved in several forms [30—32 . Here, we use the HIM to position the location 
of the final few nm gap and reduce the critical power required to burn the gap. In order 
to minimize particle contamination of the chip during the milling process, the sample 
chamber is cleaned overnight with an integrated oxygen plasma cleaner. Imaging of the 
flake is made with lower magnification in order to prevent damage before milling. 

Figure 2(a) is a HIM image of a MLG flake before milling. This image is used to 
pattern rectangles that define the areas where the flake will be milled. After defining the 
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Figure 2. (a) Helium ion microscopy (HIM) image of a MLG device before milling 
between two Au electrodes, (b) HIM image of the same device after milling. The 
resulting constriction is 100 nm wide, (c) I — 14 characteristics for the electroburning 
of the constriction in (b). (d) AFM image of the same device after electroburning. The 
inset shows the final gap across the constriction. The scale bar is 100 nm. 


milling area, the beam is focused on the edge of the flake and moved to the pattern. For 
ion beam milling an acceleration voltage of 25 kV is used with a beam current of 1 pA. 
Typical line doses are about 120 nC/cm. Two rectangular patterns, one on each side of the 
flake, become 10 nm-wide trenches after direct helium ion beam writing forming the final 
constriction in the center of the flake. In Fig. 2(b) we show the device after milling with 
the HIM. The resulting constriction is about 100 nm wide. The trench width and depth 
(^ 40 nm from AFM scan) ensure no residual connection between the milled multi-layer 
graphene edges. While the trench does extend into the oxide, we find no evidence for 
substrate damage or gate leakage after milling. 



























This is the post-peer reviewed version of the following article: 
J.O. Island et al. “Fabrication of hybrid molecular devices using multi-layer graphene break junctions”. 

J. Phys.: Condens. Matter (2014) doi: 10.1088/0953-8984/26/47/474205 

Which has been published in final form at: 
http://iopscience.iop.org/0953-8984/26/47/474205/article 






6 

4 

2 

8 

6 

4 

2 


-40 -20 0 20 40 


(b) 


V g (V) 


Figure 3. Electrical characterization after electroburning for two typical devices, (a) A 
characteristic tunneling curve is measured for the device after electroburning (red circles 
for device 1, blue circles for device 2). The Simmons model fit is plotted along with 
the tunnel current (black dashed line), (b) / — V g characteristics for the devices after 
electroburning at a bias voltage of 500 mV. 


After pre-patterning with the HIM, electroburning is performed as before but the 
constriction results in a significant reduction of the critical power required to start the 
process. Fig. 2(c) shows the / — V& characteristics of the electroburning process for the 
patterned device in Fig. 2(b). The critical power is ~ 1 mW which is more than an order 
of magnitude smaller than the narrowest as-exfoliated MLG device in Fig. 1(c). Critical 
powers for devices with wider pre-patterned constrictions are also shown in Fig. 1(c) (blue 
circles). Critical powers steadily increase to 11 mW for constriction widths up to 500 nm. 
In Fig. 2(d) we show an AFM image of the device after electroburning. The inset shows 
the gap formed across the constriction. 

We perform electrical measurements at room temperature and in vacuum to 
characterize the tunnel gap of the MLG break junctions. In Fig. 3 we show the / —14 and 
I — V g characteristics for two typical devices after electroburning. A tunneling current is 
measured through the vacuum barrier and shown in Fig. 3(a) (red and blue circles for 
devices 1 and 2, respectively). From the tunneling current, we estimate the gap distance 
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using the Simmons model for tunneling through a thin insulating him 13,33 . 


lA Ad iV) = A 


2n hd 2 


((f)- /i L )e 


-2d 


y/2 


n - (0 - hi?)e 


-2d 


y/2 m(tj>-ii R ) 


(3) 


Here, 0 is the potential barrier height, A is the cross sectional tunneling area, d is the gap 
distance, and hl = +eV/2 (/iR = —eV/2) are the chemical potentials in the left (right) 
lead. The fit using the Simmons model is plotted (black dashed line) along with the data in 
Fig. 3(a) where we have set A = 1000 nm 2 given the cross sectional area of the fabricated 
constriction (W = 100 nm, t = 10 nm). From the fit, we extract a tunnel gap distance of 
c! « 2.1 nm and a barrier height of 0 « 0.52 eV for device 1 (d « 2.1 nm and 0 ~ 0.65 
eV for device 2). The low barrier height, as compared with that of graphene (4.5 eV), 
is common for the Simmons model when applied to small junctions (e.g. Au nanogaps) 
but could also include a barrier lowering due to water or OH groups absorbed at the gap 


edge 34-36 . While the barrier height has been found to deviate from the expected value, 


scanning tunneling microscopy studies of Au nanogaps have shown a consistency between 
the extracted gap size from the Simmons fit and the actual electrode separation |37|. We 
expect that this is the case for MLG nanogaps as well but further SEM studies should be 
carried out to verify this. In Fig. 3(b) we plot the I — V g characteristics for the empty 
tunnel gap at a bias voltage of 500 mV. The current is independent of back-gate voltage. 

The high temperatures from electroburning as-exfoliated flakes are detrimental to the 
fabrication of MLG break junctions with superconducting contacts. Moreover, to achieve 


proximity effect in MLG for hybrid devices, channel lengths of < 0.5 /im are needed 38 


The combination of high flake temperature and short channel length result in melting of 
the MoRe alloy for the wider (> 3 /dm) as-exfoliated devices. This can be seen in Fig. 4(a) 
which is an AFM image of an as-exfoliated flake device after electroburning. Deformation 
of the lower contact can be seen which is near the gap shown in the inset of Fig. 4(a). By 
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Figure 4. (a) AFM image of a device without prepatterning. The high temperature 

of the electroburning process deforms the lower electrode. Inset shows the location of 
the nanogap near the deformaded electrode (scale bar 500 nm). (b) AFM image of an 
electroburned MLG break junction pre-patterned using an oxygen plasma etch. The inset 
shows the resulting constriction from etching (scale bar 100 nm). The width is ~ 80 nm. 
(c) I — 14 characteristics of the device in (b) (red circles) as well as a fit of the Simmons 
model (black dashed line), (d) I — V g characteristics of the device in (b) at a bias voltage 
of 500 mV. 


pre-patterning the MLG flake, the electroburning process can be started at much lower 
critical powers. In Fig. 4(b) we show an example of an electroburned MLG break junction 
patterned by oxygen plasma etch with MoRe contacts. The MoRe contacts are unaffected 
and the gap can be seen located at the defined constriction in the inset of Fig. 4(b). In 
Figs. 4(c) and 4(d) we show the / — V&, / — V g characteristics for the MLG electrodes with 
MoRe contacts shown in Fig. 4(b). Again, the Simmons model is fit to the tunnel current 
setting A = 800 nm 2 given the cross sectional area of the constriction (W = 80 nm, t = 10 
nm) and extract a gap distance of d ~ 1.6 nm, a barrier height of 0 1.4 eV. For the 16 

devices that resulted in tunnel gaps after electroburning, the average tunnel gap from the 
Simmons fit is 2.1 ± 0.3 nm and the average barrier height is 0.89 ± 0.3 eV. 

To demonstrate the fabrication of hybrid devices, we have performed deposition of 
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anthracene-functionalized copper curcuminoid molecules ([(Phen)CuCl(9Accm)]) on the 


MoRe MLG break junctions 39 . Fig. 5(a) shows a scheme of the molecule. Contact 
to the molecule is mediated through the interaction of two 7r-systems; the 7r-conjugated 
system of the anthracene groups on the molecule and that of the MLG break junction. For 
deposition, a powder of synthesized molecules is dissolved in dichloromethane (CH 2 CI 2 ) at 
a concentration of 0.1 mM. The electroburned junctions are kept in the solution for 30 mins 


which results in a sub-monolayer of molecules on the surface of the Si/Si02 substrate 13 


Of the remaining sixteen tunnel junctions after electroburning, one junction showed a 
significant change in the / — Vj, and I — V g characteristics after deposition. The inset of 
Fig. 5(b) shows the tunneling characteristics of the junction (low bias resistance of 14 
GO) and the inset of Fig. 5(c) shows that the current is virtually independent of gate 
voltage before deposition. After deposition, the / — 14 characteristics become linear and 
the low bias resistance is 9 MO. In addition, the current shows a strong dependence on 
gate voltage, changing by two orders of magnitude over 80 Volts on the back-gate. To 
ascertain the presence of a single molecule or a few molecules in the junction, further 
transport measurements can be performed at low temperature. 

In conclusion, we have studied the yield and stability of MLG break junctions for 
use in hybrid molecular devices. Using HIM milling or an oxygen plasma etch to localize 
a constriction in the MLG flake, we control the position of a few nm gap created by 
electroburning. The critical power needed to start the electroburning process can be 
controlled by the width of the constriction. For lower powers (~ 1 mW), we fabricate 
MLG break junctions with MoRe contacts. Using these hybrid devices, we aim to study 
superconducting effects in single molecule transport. 

The authors would like to thank Dr. Emile van Veldhoven (TNO) for his help 
with the HIM and Dr. Andres Castcllanos-Gomez for helpful discussions and manuscript 
preparation. N.A.-A. would like to thank the Spanish MCI (CTQ2012-32247). This work 
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Figure 5. MoRe MLG break junction before and after molecule deposition, (a) 
Ortep representation of the copper curcuminoid molecule [(Phen)CuCl(9Accm)]. Phen 
stands for phenantroline and 9Accm corresponds to the curcuminoid ligand that contains 
anthracene groups. In the picture the copper atom is in light blue, the carbon atoms 
are grey, the oxygen atoms are red, the N atoms are blue and the Cl atom is green. 
Hydrogen atoms are not shown for simplification, (b) I —14 characteristics before (inset) 
and after molecule deposition, (c) / — V g characteristics before (inset) and after molecule 
deposition. The / — Vj, and I — V g characteristics before deposition are also plotted in 
the main panel of (b) and (c) (blue lines) for comparison. 
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